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ABSTRACT 



Wc examine the physical processes regarding the radiatively driven mass accretion onto 
galactic nuclei due to the intensive radiation from circumnuclear starbursts. The radiation 
from a starburst not only contracts an inner gas disk by the radiation flux force, but 
also extracts angular momenta due to the relativistic radiation drag, thereby inducing 
an avalanche of the surface layer of disk. To analyze such a mechanism, the radiation- 
hydrodynamical equations are solved with including the effects of radiation drag force as 
well as radiation flux force. As a result, it is found that the mass accretion rate due to 
the radiative avalanche is given by M(r) = ri{L^. / c^){r / RY{A.R/ R){1 — e~^) at radius r, 
where the efficiency 77 ranges from 0.2 up to 1, and R are respectively the bolometric 
luminosity and the radius of a starburst ring, Ai? is the extension of the emission regions, 
and r is the face-on optical depth of a disk. In an optically thick regime, the rate depends 
upon neither the optical depth nor the surface mass density distribution of the disk. The 
present radiatively driven mass accretion may provide a physical mechanism which enables 
mass accretion from 100 pc scales down to ~pc, and it can eventually link to an advection- 
dominated viscous accretion onto a massive black hole. The radiation-hydrodynamical and 
self-gravitational instabilities of the disk are briefly discussed. In particular, the radiative 
acceleration possibly builds up a dusty wall, which shades the nucleus in edge-on views. It 
provides another version for the formation of an obscuring torus. 

Key words: accretion, accretion disks — hlack hole physics — galaxies: active — galaxies: 
nuclei — galaxies: Seyfert — galaxies: starburst — hydrodynamics — quasars: general — 
radiative transfer 
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1 Introduction 



In a standard picture of active galactic nuclei (AGNs), a viscous accretion disk of <pc 
surrounding a putative supermassive black hole is thought to be the central engine for the 
enormous energy output. The mass supply onto such tiny regions is an issue of growing 
interest. The accretion driven by a non-axisymmetric wave or a bar is a possible mechanism 
(Shore & White 1982; Shlosman, Prank, & Begelman 1989), but it is found so far that 
the bar-driven accretion is just effective in the galactic nuclear regions beyond ~100 pc 
(Wada & Habe 1995). Hence, another mechanism that enables the mass accretion from 
100 pc down to pc scales is required. Recently, Umemura, Fukue, & Mincshige (1997) have 
proposed a possible physical mechanism in the relevant scales, that is, the mass accretion 
driven by radiation from circumnuclear starbursts. 

In recent years, there have been accumulated observational evidences for starbursts in 
circumnuclear regions of AGNs. Among them, it is remarkable that molecular gas of at least 
IO^^Mq and dust of ~ lO^M© are found towards a QSO (BR1202-0725) at redshift z = 4.69 
(Ohta et al. 1996; Omont et al. 1996). This, taking into account the Universe age of ~ IGyr 
at z = 4.69 and the inferred metallicity, suggests that the QSO inhabits a massive host 
galaxy undergoing starbursts. Also, the recent HST images of nearby luminous QSOs have 
shown that luminous QSO phenomena occur preferentially in luminous host galaxies, often 
being elhpticals (McLeod &: Rieke 1995; Bahcall et al. 1997; Hooper, Impey, & Foltz 1997). 
In addition, AGN events appear to be frequently accompanied by starbursts (Soifer et al. 
1986; Scoville et al. 1986; Heckman 1991; Rieke 1992; Scoville 1992; Fillipenko 1992). The 
IRAS survey has revealed that at bolometric luminosities greater than IO^^Lq, the space 
density of starburst galaxies are just comparable to that of QSOs in the same luminosity 
range (Soifer et al. 1986). Millimeter observations of the IR galaxies indicate that all are 
exceedingly rich in molecular gas with up to 2 x IO^^Mq in their central regions (Scoville et 
al. 1986). The correlations between X-ray and CO luminosities as well as between CO and 
far-infrared luminosities in Seyfert galaxies also suggest that more powerful AGNs inhabit 
more active star- forming galaxies (Yamada 1994). Such increasing evidences suggest a 
possible physical connection between AGN events and starburst activities (Norman & 
Scoville 1988; Rowan-Robinson & Crawford 1989; Terlevich 1992; Scoville 1992; Terlevich 
1992; Perry 1992; Taniguchi 1992; Rowan- Robinson 1995). 

As for the structures of starburst regions, recent high-resolution observations of Seyfert 
nuclei, including those from HST, have revealed circumnuclear starburst regions with radial 
extension of a few tens of pc up to kpc, which often exhibit ring-like features (Wilson et al. 
1991; Forbes et al. 1994; Mauder et al. 1994; Buta, Purcell, & Crocker 1995; Barth et al. 
1995; Maoz et al. 1996; Storchi-Bergman, Wilson & Baldwin 1996; Leitherer et al. 1996). 
The rings frequently contain a number of extremely compact star clusters of < 5 pc, which 
are composed of hundreds of 0-type stars. Recently, an unexpectedly large (~ 100 pc) 
disk of gas and dust is discovered by the HST around a bright nucleus in an active galaxy 
NGC 4261 (Jaffe et al 1993; Ferrarese, Ford, & Jaffe 1996). This disk is surrounded by a 
fat torus composed of early-type stars. 
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Umemura, Fukue, & Mineshige (1997, hereafter referred to as Paper I) have considered 
the radiative effects on nuclear disks due to circumnuclear starbursts, and proposed the 
possibihty of a radiative avalanche, that is, the mass accretion driven by the relativistic 
radiation drag. In paper I, the effects of radiation drag are analyzed solely by solving 
the azimuthal equation of motion. In this paper, we explore extensively the radiative 
avalanche with including radial motion driven by radiation flux force. Also, the radiation- 
hydrodynamical instability of the nuclear region is analyzed with including vertical motion 
as well, so that the radiative building-up of a dusty obscuring wall is predicted. Some 
related topics such as the supernova effects, the generation of seed magnetic fields, and a 
possible link to an advection-dominated viscous flow are discussed. In §2, the basic equa- 
tions which govern the radiation-hydrodynamical disk evolution are provided. In §3, the 
timescales for a radiative avalanche, the disk evolution, and the resultant mass accretion 
rate are presented. In §4, the radiation-hydrodynamical/gravitational instabilities are an- 
alyzed. In §5, we discuss the effects of stellar evolution, the generation of magnetic flelds, 
and a link between a radiative avalanche and an advection-dominated accretion onto a 
black hole. §5 is devoted to the conclusions. 



2 Basic Equations 

2.1 Radiation Fields by Starbursts 

Here, we suppose a geometrically thin ring of starburst regions for simplicity, so that 
the assumption allows us to explore the radiative effects analytically. The effects of the 
extension of the starburst regions are briefly discussed. Under the assumption, the radiation 
flelds by the starburst ring with the radius R and the bolometric luminosity are given 
as follows (see also Paper I): The radiation energy density at r(< i?) is given by 

47rc(i?2 - r2) ■ ^ ^ 

The radiation flux in radial directions is given by 

L, 2x{x^ + l)Gi - 3xG2 



87ri?2 (^2 + 1)5/2 ' (2) 

where x = r/R, and Gi arc G2 are respectively Gi = G[3/4, 5/4, 1, 4x2/(a;2 -|- 1)2], and 
G2 = ^[5/4, 7/4,2, 4a;2/(a;2 + 1)2] with G being the Gaussian hypergeometric function. 
This solution is well approximated by 

,3 



47ri?2 



(3) 



2R ttR^{R - r) 

[This is slightly different from equation (2) in Paper I, but this is a more accurate approx- 
imation.] The rotation of the ring generates the toroidal flux at r to be 

= ^ (4) 
2R ' ^ ' 
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where V is the rotation velocity of the ring. The radiation stress tensors are given as 



2.2 Radiation Hydrodynamical Equations 

As a relativistic result of radiative absorption and subsequent re-emission, radiation fields 
exert drag force on moving matter in resistance to its velocity. This radiation drag extracts 
angular momentum from a rotating disk, thereby allowing the gas to accrete onto the 
center. This is known to be the Poynting- Robert son effect in the solar system (Poynting 
1903; Robertson 1937). The radiation drag usually makes effects of 0{v/c), but it could 
provide a key mechanism for momentum/angular momentum transfer in environments of 
intensive radiation fields (Loeb, 1993; Umemura, Loeb & Turner 1993; Fukue & Umemura 
1994; Umemura & Fukue 1994; Tsuribe, Fukue, & Umemura 1994; Tsuribe, Umemura 
& Fukue 1995; Fukue & Umemura 1995; Umemura, Fukue & Mineshige 1997; Tsuribe & 
Umemura 1997). 

Specific radiation forces which are exerted on moving fluid elements with velocity v are 
given as 

f^ = ^ _ y^E - VrP'-' - v^P'P') (6) 
c 

and 

= - v^E - v^P^'^ - VrP'n, (7) 

(Mihalas & Mihalas 1984) respectively in radial and azimuthal directions. Here x is the 
mass extinction coefficient, which is given by x = (k + nf.o'T)/p, where k is the absorption 
coefficient, Ue is the electron number density, ut is the Thomson scattering cross section, 
and p is the mass density. Coupled with the radiative quantities presented in the previous 
section, the last term of each component of specific forces turns out to be 0{vV/c^), and 
thus is neglected in the present analyses. Using equations (^) and @, we have the radiation 
hydrodynamical equations to 0{v/c) as 



dv^_V^_ Idp _^^Xfpr_^ 

dt r pdr dr c \ 2 



+ ±lF--:-EvA, (8) 



r dt 2c \R 

where p is the gas pressure, and $ is the external force potential (e.g., gravity). The 
azimuthal equation of motion (|^) implies that the radiation flux force (the first term on 
the right-hand side) makes fluid elements tend to co-rotate with the stellar ring, whereas 
the radiation drag (the second term on the right-hand side) works to extract the angular 
momenta from the elements. It is noted that both terms are 0{v/c) in this equation, in 
contrast to the radial flux force of 0(1) in equation (H). 
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3 Radiative Avalanche 



3.1 Timescales 

As shown in Paper I, the radiation drag time scale is given by 

First, wc consider an ionized disk. The total number per unit time of ionizing photons 
from starburst regions is roughly estimated as N\jy ~ fuYL^/hvL ~ 5.4 x 10^^s~^(L*/3 x 
IO^^Lq), with Vl being the Lyman limit frequency, and /uv is a factor which repre- 
sents the contribution to UV light by starbursts (assumed to be /uv ~ 0.1 here), while 
the total recombination rate to excited levels of hydrogen is iVrec = Ueofi^'Mg/mp ~ 
1 X 10^^s~^(Mg/10^Mo)(-Rp/100pc)~^, where Mg and Rg are the mass and the radius of 
the gas disk and a^^^ is the recombination rate coefficient to all excited levels of hydro- 
gen. The optical depth of ionized gas is roughly estimated as r ~ MgaT/T^RgiTip, which is 
numerically r ~ O.26(M(,/lO^M0)(i?g/lOOpc)~^. Thus, photo-ionization could be possible 
for an optically thin disk or an optically thin surface layer in an optically thick disk. For 
perfectly-ionized gas, the drag force is equivalent to the Compton drag, that is, x = d'T/^p, 
where ctt is the Thomson-scattering cross-section and rup is the proton mass. Then, the 
radiation drag timescale is estimated as 

(11) 

This is fairly long compared to the duration of AGN phase or starburst phase inferred from 
observations (~ lO^yr). 

In practice, a large amount of dust is detected in luminous IR galaxies (e.g. Scoville 
et al. 1991). From theoretical points of view, the timescale of photodesorption of dust by 
UV is longer than the Hubble time for refractory grains like graphites or silicates, and also 
the thermal sputtering rate is completely negligible for such grains at the temperature of 
~ lO^K (Burke & Silk 1974; Draine & Salpeter 1979). If a disk contains dust grains, then 
the accretion timescale could be much shorter because of its large opacity. The optical 
depth is estimated to be r ~ 5 x 10^ (M^/IO^Mq) (i^^/lOOpc)"^ (/^^/lO-^) , where f^g is 
the dust-to-gas mass ratio (Mathis, Rumpl, & Nordsiek 1977). If the relaxation between 
gas and dust is sufficiently rapid, then the dust accretes with accompanying gas. The 
relaxation time between gas and dust is estimated by 




STTC^my^pC^ZjlnA' 



where is the mass of dust grain, Up is the number density of protons, is the dust charge 
in units of proton charge, and In A is the Coulomb logarithm (Spitzer 1962). Tentatively 
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assuming the thickness of gas disk to be one tenth of the radius, we evaluate Up = Ue — 
1.3 X lO^cm-^(Mg/lO^°M0)(i?g/lOOpc)~^. Then, for spherical dust grains, 

W ^ 3.1 X f^V'Y^l (^] (^)\ (13) 



aO^cm-V VlO^K/ VO.l/im/ \g cm-3y V-40. 

where is the grain radius and ps is the density of solid material within the grain (e.g., 
Spitzer 1978). Hence, dust is tightly coupled with gas in a typical rotation timescale 
of ~ 10^ yr, and therefore the disk can be regarded to be composed of fluid with mass 
extinction as 

X = , (14) 

Pg+ Pd 

where n^, Cd, and pd are the number density, cross-section, and mass density of dust, 
respectively, and Pg is the mass density of gas. If taking the dust-to-gas mass ratio of 
~ 10"^ that is inferred in the Solar neighborhood, then the opacity ratio is 

ridCTd _ 1 9 X 10^ ^ «d ^ V \ ^ ( fdg \ 



UeCTT yO.l/imy ^g cm \10 J 

Thus, we find UdCTd ^ necrr and fdg <^ 1 in probable situations. Then, we can evaluate 

X - ndad/pg, so that 

«^ - ' (w^y {li^y (hT (o^) (i^) ■ (^^' 

Hence, for intensive starbursts, the drag timescale could be shorter than the duration of 
AGN phase or starburst phase. 



3.2 Disk Evolution 

Suppose that an initially quiescent gas disk is abruptly ignited by a starburst at the radius 
R with luminosity L^. Since the radiation drag efficiency decreases exponentially with the 
optical depth (Tsuribe & Umemura 1997), just an optically thin surface layer avalanches 
toward a nucleus. In equation (||), the radial radiation flux force, / c, is obviously greater 
by 0{c/v) than the drag force. So, the radiation drag makes a negligible contribution in 
the radial equation of motion. On the other hand, in azimuthal directions (see equation 
P]), both the radiation flux force and the radiation drag are 0{v/c), and therefore they 
work simultaneously. Thus, the disk evolution is expected to proceed in two steps: First, 
the optically thin surface layer (or the whole of an optically thin disk) is contracted by 
the radial radiation flux force and is quickly settled in quasi-rotation balance. Then, the 
layer sheds angular momentum due to the radiation drag, and consequently undergoes an 
avalanche. 

With the radial equation of motion (||), disregarding pressure gradient force as well as 
the drag force, the condition of rotation balance is given by 



(17) 



Here, we assume a gravitational potential such that 

dr ~ \rJ ' 

The rotation velocity of the stellar ring is given hj V = (GMdyn/ RY^"^, where Mdyn is 
the dynamical mass contained within R. If is negligibly small, then n = —1/2 simply 
corresponds to Keplerian rotation, n = to fiat rotation, and n = 1 to rigid rotation. 
Equation (|16D is rewritten with ([TtD into a dimensionless form as 

„2 



1 



X 



X] 



with using x = r/i?, t>^ = v^/V, and V = L^/L^, where is the Eddington luminosity 
corresponding to Mdyn, 

_ AncGMdyn 



X 

For a dusty disk, the Eddington luminosity is numerically given as 

/ Mdyn ' ' " 



(19) 



Le - 1.68 X 10 Lo j j ^^^^^^ ^^_3^ 

The azimuthal equation of motion (^ can be also rewritten in a dimensionless form as 

d{xv^) _ x{x - v^) 



(20) 



dt 



:i-x2) ' 



(21) 



where t = t/t-y. 

To begin with, we see the first stage, i.e., the radial contraction phase. As shown 
in equation (0), the radial compression by the flux force is expected to be strongest 
on the outer edge of disk. Because the radial flux does not extract angular momentum, 
the angular momentum conservation determines the radius of rotation equilibrium. Since 

= 1 at the outer edge of the disk before a starburst, the equilibrium radius after a 
starburst is determined by the condition of v^p = 1/x. Combined with equation (|T^), 
the resultant equilibrium radii Veq for the material that is originally located at the outer 
edge are shown versus F for three types of rotation potentials in Figure 1. The Veq is a 
slowly decreasing function with F because of angular momentum barrier. The asymptotic 
behavior for F < 1 is r^jR = 1 - [F/27r(n + 1)]^/^ while r^jR = (2/F)^/4 for F > 1. 
Equation ([T8|) shows that the radial radiation flux enhances the rotation velocity, so that 
a shear flow may be excited on the boundary between the optically thin surface layer and 
the optically thick deep stratum. The fractional difference of rotation velocity, which is 
defined by Av^/v^ = {v^ — x")/x", is also shown at r^q in Figure 1. Av^^/v^ is found to 
be weakly dependent upon F, and we have asymptotic solutions as 

(n + l)F^'/' 



Aw. 



27T 



■■F < 1 



■■F > 1. 



(22) 
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If r = 0(1) (which is indeed the case for active galaxies), Af^/f<^ = 0.2 — 0.6 in the range 
of n = —1/2 to 1. This shear flow may generate turbulence due to the Kelvin-Helmholtz 
instabilities, or the viscosity could thermalize the shear flow energy. 

Next, we consider the avalanche phase. If the radiation is very weak (F ^ 1), equation 
(^) coupled with (|T8|) is reduced to a differential equation as 

x'^^= (^-^"'') f23) 
di (n + l)a;"-2(l -x2)' ^ ' 

This equation can be analytically solved for some values of n. For instance, the Keplerian 
rotation potential (n = —1/2) leads to the solution 

i=-V^+^tan-'(^-^^]~l\og(- - 0.154. (24) 

If we take a Mestel disk (n = 0), then the solution is 

i = 1 — X — logx. (25) 

The rigid- rotation potential results in x=const. In other words, no accretion is driven in 
a rigid-rotation potential under the circumstances of very weak radial radiation flux (see 
also Paper I). If F = 0(1), however, the effects of radial flux force is important not only in 
radial contraction but also in the angular momentum transfer. When the radial flux force 
is included, the accretion is governed by 

x' = .^^"'rJr^ , (26) 



[l-x^)h{x)' 



from equations (0) and (pTj), where 



h{x) = {n + l)x^"-^ + Fx 



1 + 



X (6 — 5x) 
27r(l -x)2 



(27) 



As shown below soon, this equation allows the mass accretion even in a rigid-rotation 
potential. 

3.3 Mass Accretion Rate 

If we suppose a geometrically thin, optically thick disk, then the upper surface is irradiated 
by radiation from the upper half of starburst regions and the lower surface is irradiated by 
the lower half. Taking this into account, we can estimate the mass accretion rate as 



M = —J 2nrVrpdz 



-(1 + e-^/^) r'^^e-dn, (28) 
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where z is the vertical coordinate, dT^ = xpdz, and r is the total face-on optical depth. 
This can be rewritten in a dimensionless form using the same notations as the previous 
subsection, 



rh^ = + J xx'e-^'dn = --xx'{l - e"^), (29) 

where rh^ = M / {L^/ c^) and x' is given by (26). In an optically thick limit, 



3 . , . 

= —-XX . (30) 

Hence, the mass accretion rate is independent of the extinction coefficient as well as the 
density distribution in an optically thick disk. In comparison, in an optically thin regime, 

= —3xx't/4. This implies that the mass accretion rate is smaller in an optically thin 
disk than in an optically thick disk. This is because in optically thin media a number of 
photons penetrate the gas disk without collision with matter, and hence not all the photons 
from the starburst regions are available for radiation drag. 

The fact that the radial flux raises the rotational velocity brings two effects on the mass 
accretion rate which mutually compete: One is a positive effect that the radiation drag 
force becomes more efficient due to enhanced f<^, and the other is negative effect that the 
enhanced leads to the increase of specific angular momentum at fixed r, so that the 
accretion is decelerated. In Figure 2, for a wide variety of F = L^^/L^, the resultant mass 
accretion rate is shown as a function of radii in an optically thick regime. In this figure, a 
thin solid curve shows the estimation in the case of the negligible effects of the radial fiux, 
i.e., 

31 -n 2 

= X (-31) 



(Paper I). The thick curves are the present results, where filled circles represent the equi- 
librium radii of outer-edge shown in Figure 1. As shown in Figure 2a, for a Keplerian disk, 
the negative effect overwhelms the positive one around the outer-edge of the disk. On the 
other hand, the positive effect is predominant for a Mestel-like disk {n = 0) as shown in 
Figure 2b. In the vicinity of the ring (x — 1), we have an asymptotic solution as = 3/4, 
being independent of F as well as n. In the opposite limit, i.e., x — > 0, the solutions are not 
dependent on F for —1 < n < 1, and then the mass accretion rate is basically determined 
by (iH). 

It is worth stressing that the radial fiux can allow the accretion even if the disk is 
governed by a rigid- rotation potential (n = 1). As shown in Figure 2c, the induced mass 
accretion is not sensitive to F, and the maximal rate is M ~ 0.2L^/c^ in the range of 
F = 10^'^ to 10^. At small radii, there is an asymptotic solution as 



3 2 
= -X 



^+2 



(32) 

This approaches a limiting value 3x^/8 for F oo, which coincides with (^) for n = 0. 
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In the long run, the mass accretion rate has turned out to be almost insensitive to the 
rotation law and given by 

where t] ~ 0.2 — 1. Such radial dependence as implies that an avalanche occurs in 
a homologous manner. The radial dependence slightly deviates from near the outer 
edge of the disk. Especially, it is faster than in a rigid-rotation potential, so that some 
amount of mass would accumulate around r^q. 

Although we have assumed a ring-like configuration for starburst regions, the regions 
may be extended over a wide range of radii and heights. However, equation (|I]) shows that 
the radiation from the inner edge of the star-forming regions would be most contributive 
as long as the surface brightness falls faster than J]l{R) oc (1 — r^/i?^), which is basically 
constant for r. In other words, the regions of peak surface brightness would be always 
most responsible for the radiative avalanche. Also, in an optically thick regime, the effects 
of viewing-angles from starburst regions to a disk might reduce the mass accretion rate, 
depending upon the configuration of starburst regions (Taniguchi 1997). If we suppose a 
torus of emission regions with the half thickness of AR, then the rate should be multiplied 
roughly by AR/R due to a geometrical dilution effect (Ohsuga et al. 1997). 



4 Stability of Disk 

4.1 Radiation- Hydrodynamical Instability 

If we suppose a geometrically thin, optically thick disk, the upper surface is irradiated 
downward, while the lower surface is irradiated upward. Hence, the disk is rather confined 
by radiation force. But, the outer edge of an optically thick disk (of course, an optically 
thin disk as well) could be vertically unstable, especially if the starburst luminosity is 
super-Eddington. As far as a Kepler disk or a Mestel disk is concerned, the gravity is 
always dominant at small radii [see (^]. Thus, the inner disk remains stable, whereas the 
outer edge could be blown away by radiative acceleration. The radius inside which the 
disk is stable is given by the root of an equation as 



1 X 



2 



2 7r(l-x) 



(34) 



The solutions are shown in Figure 3 for n = —1/2, 0, and 1. We find that the radius 
is sensitively dependent upon F if r>l. To see the fate of the blown-out gas in three- 
dimensional space, we compare the radiation flux force with the gravitational force in Figure 
4, illustratively for n = —1/2. On the hatched sides for each value of F, the radiative force 
overwhelms the gravity in the vertical directions. The dotted lines are corresponding to 
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the unstable branches below which gas is accelerated downward by the gravity, while above 
which gas blows upward due to radiative force. The solid hnes are stable branches. When 
r > 2, no stable branch exists and therefore the blown-out gas may result in a superwind 
(a "blizzard" phase). [If we assume the initial mass function which is deficient in low-mass 
stars as observed in starburst regions (see also §4.1), F > 2 is corresponding to early stages 
of < 5 X 10'^ yr after the coeval star-formation. ] If the gravity by the bulge stars significantly 
decelerates the wind, the spurted-out gas may fall onto the inner disk with losing angular 
momentum by a fraction of ~ R/Vt^ (cf. Meyer & Meyer-Hofmeister 1994). When F < 2, 
the stable branch composes a wall surrounding the starburst ring. In particular, the wall 
becomes of torus shape if F < 1. Taking account of the fact that the radiation force becomes 
ineffective for optically-thick gas, the optical depth of the formed wall should be of 0(1). 
The radiation emitted originally in UV is absorbed by this wall and re-emitted possibly in 
IR wavelength. However, the radiative avalanche is almost independent of the spectrum of 
the radiation, but just depends upon the bolometric luminosity. Hence, the emission from 
the wall could induce an indirect radiative avalanche. Also, such a radiatively built-up 
obscuring wall may diminish the geometrical dilution, although the direct avalanche may 
be strongly subject to the geometrical dilution if the thickness of the starburst regions is 
as small as lOpc (Taniguchi 1997). 

The radiatively built-up obscuring wall may be a different version of the formation of 
an obscuring torus from those proposed so far (e.g. Krolik & Bcgclman 1988; Yi, Field, & 
Blackman 1994). In the present mechanism, the angular extension of the obscuring wall 
is a decreasing function of time according as the starburst luminosity fades out. Thus, 
the classification of AGN types by anisotropic viewing could be comprehended in the 
contexts of the evolution of circumnuclear starbursts. The further detailed confrontation 
with observations will be made elsewhere. 

Finally, even an inner gravity- dominant disk could be unstable to warping through a 
similar mechanism proposed by Pringle (1996, see also Maloney et al. 1996), although the 
radiative acceleration is reduced by a factor of ~ 1/r. If a radiative backreaction torque 
can drive a warp, the illumination to the warped disk is augmented, thereby increasing the 
accretion. 

4.2 Gravitational Instability 

We consider three components of gravity sources, that is, a central massive black hole, 
a massive disk surrounding the black hole, and the galactic bulge. The surface density 
distribution of the disk is assumed to be Mestel-like, and the bulge is to be uniform. To 
see the stability, we evaluate the Toomre's Q to find 

Q = 1^ + 2/, + 4/buigeX^ j ' , (35) 

where V = {GMf^yn/ RY^"^ is the virial velocity for the total dynamical mass Mdyn within 
i?, and /bh, fg and /bulge are respectively the mass fraction of the black hole, the disk, and 



12 



the bulge to the whole system. Here, we have assumed that the disk radius is comparable 
to R. If the gravity of the black hole or the bulge is predominant, then Q > 1 as readily 
seen, so that the disk can be stabihzed. However, if the coohng is effective to give Cs <S V, 
then the disk can be gravitationally unstable to a high degree. In some cases, the presence 
of an appreciable velocity dispersion could help to weaken the instability of the disk. For 
instance, the nuclear disk in Arp 220 is stabilized in outer regions of >400pc due to a 
large velocity dispersion of ~ 90km s~^, while the inner region of <400pc is marginally 
unstable (Scoville, Yun, & Bryant 1997). Anyhow, if the disk is unstable in some regions, 
it can fragment into lumps, presumably forming stars. Since the radiation drag works just 
as a surface effect for optically thick lumps, the effects of the drag force would be neghgible 
for them. Nonetheless, if at least a portion of 1/r is left in the form of diffuse gas (it is 
not unlikely because r > 10^), it may be subject to effective drag. Even if all the matter 
was processed into stars, the stellar mass loss would supply diffuse gas. Also, the newly 
formed stars may exert radiative force on the innermost gas disk. 

5 Discussion 

5.1 Sequential Starbursts 

Starbursts might not be coeval, but could be sequential inward. There are a lot of pos- 
sibilities for the sequential starbursts. The radial radiation flux can compress the outer 
edge of gas disk, where star formation may be triggered. If the UV radiation produces the 
D-type ionization front at the outer edge, massive stars could form there (Elmegreen & 
Lada 1977). If sequential explosions of supernovae occur in an OB association embedded 
in a plane-stratified gas disk, they build up an elongated hot cavity (superbubble) , and con- 
sequently the shock-heated gas spurts out from the disk (Shapiro & Field 1976; Tomisaka 

6 Ikcuchi 1986; Norman & Ikeuchi 1989). The superbubble may also compress the outer 
edge of the disk by shocks, and new stars may be born there. Furthermore, the successive 
interaction of starburst winds with the inner gas disk might be possible through a similar 
mechanism to that considered by Umemura & Ikeuchi (1987). Anyhow, if an compressed 
gas ring is ignited by star-formation, it may induce a further radiative avalanche on the 
inner gas disk. 

5.2 Effects of Stellar Evolution 

Recently, it is reported that in starburst regions the stellar initial mass function (IMF) 

is deficient in low- mass stars, with the cutoff of about mi = 2Mq, and the upper mass 
limit is inferred to be around mu = 40Mq (Doyon, Puxley, & Joseph 1992; Chariot ct 
al. 1993; Doane & Mathews 1993; Hill et al. 1994; Blandle et al. 1996). Here, we 
consider the effects of stellar evolution. We assume a Salpeter-type IMF for a mass range 
of [mi, mj as 0(m*) = (po{m^ / Mq)^^^~^^\ where s — 1.35 and m* is mass of a star. We 
employ simple relations of the stellar mass to luminosity as {U/Lq) — {m^/MQy with 
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q = 3.7, and the mass to age as = 1.1 x lO^°yr(m*/M0)-'^ with u = 2.7 (Lang 1974). 
Then, we obtain the total stellar mass and luminosity respectively as M* = 1.460o-^© and 
= 2.20o(82if — 1)Lq, where tj is the elapsed time in units of lO^r after the ignition 
of a starburst. These give a large himinosity-to-mass ratio in early stages of evolution 
as L*/M* ~ 130tf "'^'^(Lo/Mq), Also, if we assume that stars of > AMq are destined to 
undergo supernova explosions and release the energy radiatively with the efficiency of csn to 
the rest mass energy, the temporally averaged supernova luminosity is Lg^ = esN-^SNC^ = 
2000^7 "'^^-f'©; where MgN(t) is the total mass of SNs which exploded until t, and esN = 10~^ 
is assumed. Thus, the ratio of Lsn to L* is Lsn/-^* = 0.12/(1 — O.Oliy^'^). Hence, the 
luminosity is always dominated by stellar luminosity until the supernova explosions are 
terminated at tr — 26. 



5.3 Generation of Seed Magnetic Fields 

Although the Compton drag is not predominant in a dusty disk, it is significant in the 
sense that magnetic fields could be generated, because the Compton drag force is exerted 
virtually on free electrons (Mishustin & Ruzmaikin 1971; Zeldovich, Ruzmaikin & Sokoloff 
1983). The equation of motion for electrons is 

- {e/meC){cE + VeX B)- U^jVe + Vep{Vp " v^) = 0, (36) 

where rrig is the electron mass, v^. and Vp are the electron and proton velocities, E and 
B are the induced electric and magnetic fields, v^-i = 'nip/rnj:^ ~ SaTL^/STmieC^R^, and 
Uep — 4.8neT'~^''^ In A (s~^) with In A being the Coulomb logarithm (~ 30). Since is 
estimated as ~ 1.2 x 10~^^s~^(L*/3 x lO^^L0)(i?g/lOOpc)~^, u^y -C t'ep, and therefore 
Vg is approximately equal to the fiuid velocity {Vg — v). Then, Maxwell equations give 

— = X - — f V^B. 37) 

Here, one may ignore the second term of the right-hand side for fields of galactic scale, 
because it arises from the Coulomb interactions (Zeldovich, Ruzmaikin & Sokoloff 1983). 
Thus, B is evaluated by 

vn c f 

B = J i/e^V X vdt. (38) 
If the gas disk is in Kepler rotation, then |V x ~ v^/r, so that 



V 1/2 / ^ \ -3/2 
'0/ 




(39) 

Remembering that the radial radiation fiux produces shear on the disk surface, these seed 
magnetic fields could be quickly enhanced due to the combination of differential rotation, 
turbulence, and magneto-rotational shear instabilities (Balbus & Hawley 1991; Hawley, 
Gammie & Balbus 1995; Matsumoto & Tajima 1995; Matsumoto et al. 1996), which may 
be relevant to inner viscous flow. 
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5.4 Link to an Advection-Dominated Flow 



Recently, the solution of an advection-dominated viscous flow has been developed in the 
prescription of a-viscosity (Abramowicz et al. 1995; Chen et al. 1995; Narayan & Yi 1994, 
1995; Narayan 1996; Narayan, Kato, & Honma 1997; Chen, Abramowicz, & Lasota 1997). 
In the present radiation-hydrodynamical model, the strong shear flow which is comparable 
to rotation velocity is induced on the disk surface, and hence a hot corona with roughly 
virial temperature may be generated. This is a favorable situation where an advection- 
dominated accretion is raised. The advection-dominated solution exists when the mass ac- 
cretion rate is relatively low as M < Merit- If a fiducial efficiency of 10% is assumed (Frank, 
King, & Raine 1992), then the critical rate is given as Merit — 3a^LE,Th('^/100rg)~^/^/c^ 
at r > lOOr^ and constant at r < lOOr^, where a is the viscosity parameter (a < 1, 
Shakura & Sunyaev 1973), Vg is the Schwarzshild radius of the central black hole and L^^Th 
is the Eddington luminosity corresponding to the Thomson scattering. The energy conver- 
sion efficiency is e = Ladv/Mc^ ~ 0.1(M/Mcrit), where 

Ladv is the energy release from an 
advection-dominated accretion onto a putative massive black hole. When equating Merit 
with (|33D , we find a transition radius where the mass accretion by the radiative avalanche 
links to an advection-dominated flow. The transition radius is given by 

where Mbh is the mass of a central black hole. [It is noted that more precisely we should 
employ the total dynamical mass in relevant scales instead of Mbh (Mineshige & Umemura 
1996, 1997), although M^jn ~ Mbh if a;tr ^ 1-] Thus, if i? ~ lOOpc, the radiative avalanche 
may link to an advection-dominated accretion around ~ Ipc. 



6 Conclusions 

We have investigated fundamental physical processes of radiatively driven mass accretion 
due to circumnuclear starbursts by solving the radiation-hydrodynamical equations of mo- 
tion including radiation drag force as well as radiation flux force. The present results are 
summarized as follows. 

1. The mass accretion rate via a radiative avalanche is almost insensitive to a gravi- 
tational potential, if the starburst luminosity is comparable to or greater than the 
Eddington luminosity. In particular, the radial flux force allows the disk gas to ac- 
crete even in a rigid-rotation potential. The mass-accretion rate at radius r is given 
by M(r) = ?7(L*/c^)(r/_R)^(Ai?/i?)(l — e""^), where rj ranges from 0.2 up to 1, and 
R are respectively the bolometric luminosity and the radius of a starburst ring, Ai? 
is the extension of the emission regions, and r is the face-on optical depth of a disk. 
In an optically thick regime, the rate depends upon neither the optical depth nor 
surface mass density distribution of a disk. The timescale of the radiative avalanche 
could be shorter than 10®yr for a dusty disk. 
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2. The radiation- hydrodynamical instability could prevent the disk surface from accre- 
tion especially in early stages ( < 5 x 10'' yr) of stellar evolution. In later stages, it 
would result in forming an obscuring wall of torus shape. The wall indirectly induces 
a radiative avalanche with diminishing the geometrical dilution (AR/R). Also, such 
a radiatively built-up wall shades the nucleus in edge-on views, so this is another 
version of the formation of an obscuring torus. In the present mechanism, the angu- 
lar extension of the obscuring wall is a decreasing function of time according as the 
starburst luminosity fades out. Thus, the classification of AGN types by anisotropic 
viewing could be comprehended in the contexts of the evolution of circumnuclear 
starbursts. This picture will be confronted with observations in further detail else- 
where. 

3. For an ionized surface of disk, since the Compton drag works virtually on free elec- 
trons, magnetic fields could be generated. Although the magnitude is as small as an 
order of ~ 10~^^ Gauss, it could be amphfied quickly by the combination of differen- 
tial rotation, turbulence, and magneto-rotational instabilities, which may be relevant 
to the inner viscous fiow. 

4. Since the strong shear fiow induced by radiation force would result in hot coronae 
with roughly the virial temperature, a radiative avalanche is likely to link to an 
advection-dominated accretion at a point where the mass accretion rate falls to the 
critical rate, e.g., < Ipc. The link to an advection-dominated fiow is possibly related 
to X-ray properties in ultra-luminous IRAS galaxies. 
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Figure 1: The equilibrium radius, req, of the outer edge of a disk after contraction by the 
radial radiation flux is plotted against F = L^jL^ {upper panel) , where is the bolometric 
luminosity of a starburst ring and Le is the Eddington luminosity corresponding to the 
gravitational mass within the ring radius, R. Here, three types of gravitational potentials 
are assumed, i.e., n — —1/2 (Keplerian rotation potential), n — (flat-rotation potential), 
and n = 1 (rigid- rotation potential). Also, the fractional difference between an enhanced 
rotation velocity by the radiation flux and an original rotation velocity of gas disk, Av^/v^, 
is shown at Teq {lower panel). 
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Figure 2: (a) The mass accretion rate by a radiative avalanche is shown in units of L^^jc? 
as a function of radius normahzed by the ring radius, i?, for F = L^/Lg = 10 {dot-dashed 
curve), 1 (dotted), 0.1 (dashed), and 10^^ (solid). Here, an optically thick disk is assumed 
in a gravitational potential of Keplerian rotation (n = —1/2). Filled circles represent the 
equilibrium radii of the outer-edge shown in Figure 1. A thin solid curve is the previous 
estimation with negligible effects of the radial flux: M/(L^/c^) = (3/8)(r/i?)^[(l -n)/(l + 
n)] (see Paper I). In the vicinity of the ring (x ^ 1), all solutions approach an asymptotic 
value of M = 3L*/4c^, being independent of F as well as n. 
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Figure 2: (b) Same as Figure 2a, but for n = (fiat-rotation potential). 
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Figure 2: (c) Same as Figure 2a, but for n = 1 (rigid-rotation potential). Although the 
rigid rotation (n = 1) results in no mass accretion if a radial flux is negligibly small, the 
radial flux for F > 10~^ can generate an appreciable accretion. 
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Figure 3: The radius inside which the disk is vertically stable is shown against F = I/^/Lg 
for n = —1/2, 0, and 1. Outside the radius, an optically thin disk or the outer edge of an 
optically thick disk could be blown out by the radiation flux force. 
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Figure 4: Vertical radiation flux force is compared with the vertical gravitational force for 
n — —1/2 (a Keplerian disk) in cylindrical coordinates (r, z). On the hatched sides for 
each value of F = L*/Le, the radiation force overwhelms the gravity. The dotted lines 
are corresponding to the unstable branches below which gas is accelerated downward by 
the gravity, while above which gas blows upward due to radiative force. The solid lines 
are stable branches. The attached ages are that of starburst regions corresponding to each 
value of F. 
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